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Abstract. We give a simple proof that the Neumann coefficients of surface states in 
Witten’s SFT satisfy the Hirota equations for dispersionless KP hierarchy. In a similar 
way we show that the Neumann coefficients for the three string vertex in the same theory 
obey the Hirota equations of the dispersionless Toda Lattice hierarchy. We conjecture 
that the full (dispersive) Toda Lattice hierachy and, even more attractively a two-matrix 
model, may underlie open SFT. 


1 Introduction 

A remarkable integrable structure underlies Witten’s String Field Theory (SFT), [jl], It 
manifests itself in the Neumann coefficients, which can be shown to satisfy the Hirota 
bilinear equations. This is true for the Neumann coefficients of the surface states such as 
the sliver and the butterfly states, [j2j, (J ft. j5|. I 11 this case the relevant integrable hierarchy 
is the dispersionless KP (dKP) hierarchy. But, what is more important, it is also true for 
the Neumann coefficients that define the three string vertex, a structure which is at the 
core of Witten’s SFT since it describes the basic interaction and the star product. At the 
tree level in this case the relevant integrable hierarchy is the dispersionless Toda lattice 
(dTL) hierarchy. 

Following an inspiring suggestion of f(J, in this paper we show in a straightforward 
way that the Neumann coefficients for surface states satisfy the Hirota equations for 


the dKP hierarchy and that the three string Neumann coefficients satisfy the Hirota 
bilinear equations for the dTL hierarchy. Since the dKP hierarchy is a subcase of the dTL 
hierarchy, it is evident that SFT is based (at the tree level) on the latter. 

These results are simple verifications that some identities hold for Neumann coefficients 
in SFT, but they may have important implications. For it is naturally tempting to 
conjecture that the full SFT is based on the (dispersive) TL hierarchy. And since the 
latter is the hierarchy underlying two-matrix models, one is lead to conjecture that at 
the basis of SFT there may lie a two-matrix model. So the discovery of the integrable 
structure of Witten’s string field theory opens the way to a new (old) exciting field of 
research. 

The paper is organized as follows. In section 2 we introduce the Hirota equations for 
dKP and dTL hierarchies. In section 3 we introduce new compact formulas for Neumann 
coefficients for surface states and verify that they satisfy the Hirota equations for the dKP 
hierarchy. In section 4 we introduce generating functions for the Neumann coefficients of 
the three string vertex and verify that they satisfy the bilinear Hirota equations for the 
dTL hierarchy. Section 5 is devoted to comments and future prospects. 


2 dKP and dTL Hirota equations 


The essence of integrable hierarchies is contained in the Hirota bilinear equations 0 . We 
present them here for the dKP and dTL hierarchies |8], |] . To this end let us introduce 
the flow parameters to? tk, tk with k — 1, 2,.., oo and the differential operators 


OO 


m = £ 

fc=i 


_LA 

kz k dt k ’ 


OO 


m = £ 

k =1 


1 d 
kz k dt k ' 


( 1 ) 


Next we need the r-function (free energy) F = F({t 0 ,t k ,t k }). For the dKP F depends 
only on {t k }, k > 1. Then the Hirota equations for the dKP hierarchy can be written in 
the following compact elegant form (for more details, see [0.0 and references therein): 

(*! - z 2 )e D{z 1)D(Z2)F + (z 2 - z 3 ) e D ^ D ^ F + (z 3 - Zl )e D(z3)D(zi)F = 0. (2) 


The dispersionless Toda Lattice hierarchy depends instead on all the parameters {t 0 , t k , t k }. 
The Hirota equations can still be written in compact form 


(z t - z 2 )e D ^ D ^ F = Zl e-°'° D ^ F - z 2 e- dt ° D ^ F , 

Zl z 2 (l - e -D(zi)D(z 2 )F\ = e d t0 (d t0 +D( Z1 )+D(z 2 ))F 


( 3 ) 

( 4 ) 


and their bar-counterparts. 

Our task in this paper is to identify the exponents in the above equations (derivatives of 
F) with suitable generating functions of Neumann coefficients and verify that the relevant 
Hirota equations are identically satisfied. 
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3 Neumann coefficient for surface states 


Neumann coefficients for surface states were defined in [I2|. Given a conformal map f(z) 
from the upper semidisk, they are defined by 


nL = 


1 


dz 1 r dw 1 


nm Jo 2ni z n Jo 2ni w m 


d z d w \n (f(z) - f(w )) 


(5) 


In the following we actually only need that f(z) admits a Taylor expansion around z = 0, 
with /'(0) 0. In (|5|) the integration contours are counterclockwise about the origin. 

Now we define the generating function N^(zi,z 2 ) for large z\ and z 2 , 


OO OO "I "I 

A^ 1 ,z 2 ) = ££--Ai r 

n =1 m= 1 "1 ^2 


( 6 ) 


Using (|^), we make this definition more explicit as follows 
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The series converge for \z\w\ > 1 and \z 2 C,\ > 1. The integration contours are chosen in 
such a way as to surround the logarithmic cuts. Next we insert the decomposition 


d w ln(l - —) <9 c ln(l - -E) = -E 


+ 


Z\W 


z 2 ( w C w(-Z 2 C W - z 1 1 w-z 1 1 ( -z 2 
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( 8 ) 


We see immediately that the term W gives rise to a divergent term in ([7|). For this 
reason we regularize N^(zi,z 2 ) by subtracting 


JUlJLd Jn(l 
0 2 7ri Jo 2m 


—) <9 c ln(l - -E) ln(w - C) 
ziw z 2 Q 


from it. Repeating the derivation, inserting (j^) and performing the contour integrals, we 
get 


N f (z 1 ,z 2 ) = 

= In 


dw r d( 


1 v 


0 2m Jo 2m n ~^ =l n(ziw) n m(z 2 () v 

m /(w 1 ) - /&- 1 ) 


d w d(\\i 


7H-/(0 ’ 


71 - ^2 (/(0) - f(z 2 ))(/(«! )-/( 0 )) 


(9) 


We stress that the subtracted term with ln(u> — C) does not change the value of 
the initial integral because it corresponds to the identity mapping with obviously trivial 
Neumann coefficients. 
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Now we identify N f (z 1: z 2 ) with D(zi)D(z 2 )F in eq.(|), 

D( Zl )D(z 2 )F = N f {z u z 2 ), 

and analogously for the other exponents. It is now elementary to prove that 
tically satisfied. We simply remark that N^(z\,z 2 ) (]9[) can be rewritten as 

1 1 \\ 


jv^(^ 1jS2 ) = 1.1 f 

V-l -Z 2 \( 


( 10 ) 

is iden- 

( 11 ) 


xm-fte 1 )) (m-Hzr 1 ))))' 

then plug this equation in (||) and use cyclicity. 

The fact that N^(zi,z 2 ) satisfies the Hirota equations means that the corresponding 
Neumann coefficients Nl m (^) satisfy all the identities implied by (^), which can be 
obtained by expanding the latter in powers of Zi, z 2 and Z3, and equating the corresponding 
coefficients. The first few of these identites have been written down in HTTI. |6fl and will not 


be repeated here. Instead we would like to point out that formula (|TT]) is very efficient in 
explicitly computing the relevant Neumann coefficients. A very simple example is given 
by the so-called Nothing state, ||. In this case 

z 


f(z) = fn(z) = 


( 12 ) 


1 + z 2 

Plugging in this formula one can obtain the following expression for the generating func¬ 
tion N^ n (zi,z 2 ) of the Neumann coefficients of the Nothing state 


N f "( Zl ,z 2 ) = 111 ( 1 -—) 

V ziz 2 J 


(13) 


Then, the Neumann coefficients can easily be extracted by expanding N* n (zi,z 2 ) using 


N fN (z l7 z 2 ) = -E E iwT 

Zi 


J n,mi 


Nil = 


J n,m 


(14) 


1=1 m=1 -61 -*2 V nm ' . y/n™ 

and, up to normalization, reproduce correctly the corresponding coefficients given in ||. 
As a second example, we mention the series of states characterizing by [5J 

/ z n ~^~^ d \ z 

/ fl „M=exp(-— -)z= (1 + Z „ )V „ n — 1,2,... . (15) 

The state at n = 2 corresponds to the so-called Butterfly state, while the states with other 
n are generalizations of the latter. Plugging © into eq. ( |TT| ) we obtain the corresponding 
generating functions N^ Bn (z\, z 2 ) of Neumann coefficients for these states 

(l + ziy/n- (1 + z2) 1 / n \ 


jy/sn Z2 ) — } n 


Zi - 2:2 


(16) 


f B 

The Neumann coefficients Nlm for the Butterfly state were calculated explicitly in 
and it is a simple exercise to verify that they are reproduced from eq. ( |T6| ) at n = 2. Let 

fB 

us remark that, actually, the method, developed in H to evaluate Nlm , can also be used 

£ 

to evaluate the Neumann coefficients Nlm for any n, and we have explicitly verified that 
they are indeed reproduced by our simple formula (|IG|). 
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4 Neumann coefficients for three string vertex 


The three string vertex 100 of the Open String Field Theory is given by 

l^s) = f d 26 p {1) d 26 p {2) d 2e p {3) S 2e (p {1) +p( 2 ) +P( 3 )) exp (-E) |0,p)i 23 (17) 

where 

e = E [\ £ rf'C^ + Evi(yrf" + 

r,s=l \ m,n>l n> 1 



Summation over the Lorentz indices n, zz = 0,..., 25 is understood and r] lw denotes the flat 
Lorentz metric. The operators denote the non-zero modes matter oscillators 

of the r-tli string (r = 1,2,3). p l [ r) are the center of mass momenta of the r-tli string. 
|p( r )) is annihilated by the annihilation operators a and is eigenstate of the momentum 
operator pjL with eigenvalue p^ r y 
The Neumann coefficients VI 


rrs 
' mn i 


v rs 

V 0 n 


and VII 


were explicitly computed in [ 131. For 


definiteness, we refer to the coefficients contained in Appendix A of [15] as the standard 
ones. 

Our purpose in this section is to show that these coefficients too obey the Hirota 
equations. The difference with the previous section is that these Hirota equations are 
the ones characteristic of the dispersionless Toda Lattice hierarchy. We recall that in the 
definition of V^ there is a gauge freedom, since, due to momentum conservation we are 
allowed to make the redefinition 


V rs —► V rs -i_ A s 

' / 0n ' / 0n ' ^n • 


(19) 


We will see below that the Hirota equations involve only gauge invariant combinations of 
them. 

The strategy is the same as in the previous section. We first define a more efficient 
method to calculate the Neumann coefficient, by introducing suitable generating functions. 
A basic difference with the previous section is that in this problem, instead of one, we 
have three functions involved 


We define 


1 — iz\ 
1 + iz 2 \ 5 
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Now, we insert in these definitions the formulas for N™ n and iVg* given by Jl6|, and, as in 
the previous section, we extract the corresponding generating functions. The task is easy 
for N rr , because we can limit ourselves to taking the formula of the previous section. In 
other words we have 


N rr (z 1 ,z 2 ) = In 


fr(z i ) - f r {z 2 ) 


' m _ 

z\ - Z 2 (fr( 0 ) - /r(^2 1 ))(/r(^r 1 ) “ /r( 0 )) 


(23) 


For r^swe start again from ||T 6 

N ~ =-—i / A J_ 

nm Jo 2ni z n Jo 2ni w m 


d z d w ln (f r (z) - f s (w )). 


(24) 


Next we insert this formula in fl2ip . We notice that, since implies f r ( 0) 7 ^ f s ( 0), we 
do not need to do a subtraction as in the previous section, so the method works straight 
away. The result is 


N rs (zi, z 2 ) = —In 


(frjzT 1 ) - jsi^mm - my 

{fr{ 0) - /s(- 2 _1 ))(/r(-r 1 ) - fs(0)), 


(25) 


As for Nq s (z ) we start again from the formula that can be found in |L 6 


N rs = 

iv 0 n 


dz 1 
m Jo 2ni z n 


9Jn(f,.(0) - f,(z)), 


(26) 


and plug it into eq. (| 22 |). Once again we have to distinguish two cases r = s and r J s. 
The first case needs a subtraction, the second does not. The result is 


N r 0 r (z) = In \z 

while, for the case r^s, we obtain 


fr( 0 ) ~ fr(z 

-/;( 0 ) 


(27) 


N r 0 s (z) = In 


7 r( 0 ) - f s (z 


- 1 ' 


fr( 0 ) - fs( 0 ) ) ‘ 

From the above definitions it is immediate to see that 


(28) 


N rs — N sr 

nm mn ’ 


Moreover, from the explicit form of the functions f r (EQ), one can easily verify that 


N r+hs+ 1 = N rs , IV 0 r+1 ’ s+1 = Nq S 


with r and s defined mod 3. 

1 Notice the different sign with respect to | 1 G|. 
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Now, having all the generating functions at hand, we want to show that they satisfy 
the Hirota equations (0-0) • To this end we proceed to identify the exponents of these 
equations with the previous generating functions. When zero modes are not involved 
things are straightforward. We have, for instance, 

D{ Zl )D(z 2 )F = N 11 (z u z 2 ), (29) 

D(z 1 )D(z 2 )F = N 12 (z 1 ,z 2 ). (30) 


The identification for the generating functions involving zero modes is more compli¬ 


cated and can only be by trial and error, because of the gauge freedom (pT|) mentioned 
above. After some attempts we came to the following identification: 

3„flMF=jV"W-CM. (31) 

d ta D(-z)F = NZ\z)-Nl\z). (32) 

We see that the combinations in the RHS are gauge invariant! 

Now we can proceed to verify that our generating functions satisfy the Hirota equa¬ 
tions. Substituting (|29|) in the LHS of (0) and (31) on the RHS, it is trivial to verify that 
the equations identically satisfied. Similarly, for eq. (0) we have to insert (j30[) in its LHS 
and (fTT] (10) in its RHS. However, in this case, we have an indeterminate: F 0 o := d 2 Q F. 
After some elementary algebra one finds that eq. (0) is identically satisfied provided 

mf’( o)_i6 

27' 


e F °° = 


(/ i ( 0 )-/ 2 ( 0)) 2 

’oo is identical to the standard coefficient 


(33) 


t rrr 

' ^00 * 


We see that F\ 

The Neumann coefficients that one can extract by expanding the generating function 
(^3|) and (|25| ) coincide, up to normalization, with the standard Neumann coefficients V r 
For instance, for r fy s, let us set N 12 (l/x, 1 /y) = H(x,y). Then 


rs 

nm 


v 12 = 

nm 


n!m ! 


-a:arH( o,o). 


(34) 


Here are some examples 
16 


V,\ 2 = 


li 


V 12 - -V 12 - - 
v 12 ~ V 21 ~ 


32 


y 12 = 16 

31 y/3 3 5 ’ 


V 2 \ 2 = yft- 


160 


27’ 1Z Zi 81V 3’ ' di y/3 3 5 ’ ’ 3 6 ' 

More generally, on the basis of the above equations, we have the following identifica¬ 
tions among second derivatives of F, standard Neumann coefficients and the Neumann 
coefficients N r n f m . 

d d 


F t t = 

L n L m 

F t t = 

f'n Lm 


dt n dt r 
d d 
dt n dt r 
= d_^_ 
totn _ dt 0 dt n 
d d 
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F = -y/nmV^n = nmN x n 
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F = Vf iy«n - C) = n (Nil - "Z) . 


' nmV nm = nmN nm , 


Ft0tn ~ dtn dtn 


F = 


— (V 21 — W llN ) — n (N u — N 21 ) 

r , (/On v 0n ) — IL \ 1 ' l 0n Jv 0 n)' 


(35) 

(36) 

(37) 

(38) 
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The Neumann coefficients involving zero modes and computed by expanding (|27|) and 
(p8|), do not coincide in general with the standard ones. This is no surprise since the latter 
were calculated with a definite gauge choice. As far as the coefficients are concerned 
we see once again that the Hirota equations only constrain gauge invariant combinations 
of them. 

The Hirota equations generate infinite many relations among the coefficients N^, 
and F 0 o, relations which one can derive by expanding (||-[|) in powers of zi,z 2 ,z 2 . For 
instance, from eq. © we get 


^titi — ~^t 0 t 2 — 2^ totl ^ 2 ’ (39) 

2^1*2 = g-Ptota + 2 ^i^ 2 — Q^toh) 3 , (40) 

(41) 


and, from ©, 


e F ” = F u r „ 

e Fou Uofi — 2-Ui (a. 


(42) 

(43) 


and so on. All these equations are identically satisfied. 


5 Discussion 


From the above results there can be no remaining doubt that an integrable structure 
underlies Open String Field Theory. At the tree level this is the dTL hierarchy, which 
implies in particular that the Neumann coefficients are powerfully constrained. This result 
in itself may seem at first simply formal, since it was already known how to calculate all 
these coefficients (but above we have presented new effective formulas to compute them). 
However this fact opens the way to new possibilities, which seemed to be out of reach 
until now. For, on the one hand, the dispersionless TL hierarchy is a limiting case of 

H 


the (dispersive) Toda Lattice Hierarchy (TLH) 


This hierarchy describes all genus 

0 , 


corrections to the dTL results. On the other hand, we know that two-matrix models, 
are precisely realizations of TLH, in the sense that any solution of two-matrix models 
satisfy the Hirota equations of TLH. We recall that the Hirota equations are the same for 
all two-matrix models, but any two-matrix model is characterized in addition by specific 
coupling conditions (or string equations). We do not actually know whether any solution 
of TLH can be incorporated in a definite two-matrix model. However this seems to be very 
likely, and, anyhow, we have indirect evidence that two-matrix models have to do with 
string theory. In the past specific two-matrix models were shown to describe topological 
held theories chracterizing string vacua with central charge c < 1, 


19 , or c = 1, 


It is therefore hard to refrain from conjecturing that a specific two-matrix model underlies 
open SFT. If true, this would be a powerful tool to extract exact higher genus results, 











and perhaps also exact all genus results. In fact we recall that, while the Hirota equations 
generate relations among the Neumann coefficients (but do not completely determine 
them), a two-matrix model would allow us, at least in principle, to fully determine them 
(see, for instance, pi Ell] for some examples in an admittedly very simple case). 


Needless to say all we said above rings a bell: two-matrix models are coming back 

Perhaps the parallelism 


via the Dijkgraaf-Vafa description of Yang-Mills theory, see j£2 


between the two appearances of matrix models is not accidental. 

Another line of research the results of this paper open up involves their supersym¬ 
metric extension, which is of course relevant for supersymmetric SFT. In this direction 
the relevant object is likely to be the dispersionless version of the supersymmetric KP 
hierarchy proposed in p|. Quite recently, the N = (111) supersymmetric dTL hierarchy 
was proposed in ][h|] which can also be relevant for this problem. However we do not 
have yet a full fledged version of the latter hierarchy, e.g. its bilinear form (if any), so the 
research in this field faces an additional obstacle and requires further investigations. 
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